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We explore a model in which the cosmological constant A and the density 

contrast at the time of recombination a^ec are random variables, whose range 

> 
(^ ' and a priori probabilities are determined by the laws of physics. (Such models 

CN ' arise naturally in the framework of inflationary cosmology.) Based on the 
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assumption that we are typical observers, we show that the order of magnitude 
coincidence among the three timescales: the time of galaxy formation, the 
time when the cosmological constant starts to dominate the cosmic energy 
density, and the present age of the universe, finds a natural explanation. We 
also discuss the probability distribution for CTree, and find that it is peaked near 



^ ' the observationally suggested values, for a wide class of a priori distributions. 

I. INTRODUCTION 

During the past year and a half two groups have presented (independently) strong evi- 
dence that the expansion of the universe is accelerating rather than decelerating |jl| . This 
surprising result conies from distance measurements to more than fifty supernovae Type la 
(SNe la) in the redshift range z = Q io z = 1.2. While possible ambiguities related to 
evolution, and to the nature of SNe la progenitors still exist 0, the data are consistent 
with the cosmological constant (or vacuum energy) contributing to the total energy density 
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about 70% of the critical density {Q\ ~ 0.7). 

At the same time, other methods, and measurements of the anisotropy of the cosmic 
microwave background indicate that matter alone contributes about Qm — 0.3, which when 
combined with the cosmological constant suggests a flat universe 0. 

These findings raise however an extremely intriguing question. It is difficult to under- 
stand why we happen to be living in the first and only time in cosmic history in which 
Pm ~ Pa (where pm is the matter density, and pA the vacuum energy density associated 
with the cosmological constant). That is, why 

to~U , (1) 

where to is the present time and ^a is the time at which the cosmological constant starts 
to dominate. Observers living at t ^ tA would find Qm ~ 1 (^a ~ 0), while observers at 
t ^ tA would find f^A ~ 1 {^M ~ 0). 

There is another, less frequently discussed "coincidence" , which also calls for an explana- 
tion. Observationally, the epoch of structure formation, when giant galaxies were assembled, 
is at 2; ~ 1 — 3, or t^ ~ ^o/3 — to/8- For the value of A suggested by observations, this is 
within one order of magnitude of ^a, 

tc ~ U. (2) 

It is not clear why these seemingly unrelated times should be comparable. We could have 
for example to <^ ^a- 

In the present work, we explore whether the above "coincidences" [eqs. (|I]) and @] 
could be due to anthropic selection effects. The approach that we use is one in which it is 
assumed that some of the constants of nature are actually random variables, whose range 
and a priori probabilities are nevertheless determined by the laws of physics. Under this 
assumption, some values which are allowed in principle, may be incompatible with the very 
existence of observers. Hence, such values of the constants cannot be measured. The values 
in the observable range will be measured by civilizations in different parts of the universe, and 



we can define the probability dV = P{x)dxi---dxn for variables Xa to be in the intervals dxa 
as being proportional to the number of civilizations that will measure Xa in those intervals. 
Following Ref. |Q, we shall use the "principle of mediocrity", which assumes that we are 
"typical" observers. Namely, we can expect to observe the most probable values of Xa- 

An immediate objection to this approach is that we are ignorant about the origin of 
life, let alone intelligence, and therefore the number of civilizations cannot be calculated. 
However, the approach can still be used to find the probability distribution for parameters 
which do not affect the physical processes involved in the evolution of life. The cosmological 
constant A and the amplitude of density fluctuations at horizon crossing Q are examples of 
such parameters. If the parameters Xa belong to this category, then the probability for a 
carbon-based civilization to evolve on a suitable planet is independent of Xa, and instead 
of the number of civilizations we can use the number of habitable planets or, as a rough 
approximation, the number of suitable galaxies. We can then write 

P(x)rf"x « dM, (3) 

where dM is the number of galaxies that are formed in regions where Xa take values in the 
intervals dXa- 

The problem of calculating the probability distribution dV(x) can be split into two 
parts. The number of galaxies dAf{x) in Eq. (^ is proportional to the volume of the 
comoving regions where Xa take specified values and to the density of galaxies in those 
regions. The volumes and the densities can be evaluated at any time. Their product should 
be independent of the choice of this reference time, as long as we include both galaxies that 
formed in the past and those that are going to be formed in the future. For some purposes 
it is convenient to evaluate the volumes and the densities at the time of recombination, tree- 
We can then write 

dVix) = y{x)dV.{x)- (4) 

Here, dV^:{x) = P*{x)d^X is proportional to the volume of those parts of the universe where 
Xa take values in the intervals dxa, and i^ix) is the average number of galaxies that form per 



unit volume with cosinological parameters specified by the values of Xa- dV^^x) is an a priori 
probability distribution^ which should be determined from the theory of initial conditions 
(e.g., from an inflationary model). On the other hand, the calculation of z^(x) is a standard 
astrophysical problem, unrelated to the calculation of the volume factor dV,^{x)- 

The principle of mediocrity (which is closely related to the "Copernican principle") has 
been applied to determine the likely values of the cosmological constant [^|-0] , of the density 
parameter Q [l|,i|, and of the density fluctuations at horizon crossing Q [0. A very similar 
approach was used by Carter |Tl[] , Leslie [|l^ and Gott |]l3l to estimate the expected lifetime 



of our civilization. Gott also applied it to estimate the lifetimes of various political and 
economic structures, including the journal ^^Nature" where his article was published. Related 
ideas have also been discussed by Linde et. al. ||1^ and by Albrecht |jl5| . 



Spatial variation of the "constants" can naturally arise in the framework of inflationary 
cosmology [0. The dynamics of light scalar fields during inflation are strongly influenced by 



quantum fluctuations, causing different regions of the universe to thermalize with different 
values of the fields. For example, what we perceive as a cosmological constant could be a 
potential f/(0) of some field 0(x). If this potential is very flat, so that the evolution of 
is much slower than the Hubble expansion, then observations will not distinguish between 
[/(0) and a true cosmological constant. Observers in different parts of the universe would 
then measure different values of U{(j)). Quite similarly, the potential of the inflaton field 
$ that drives inflation can depend on a slowly-varying field 0. In this case, regions of the 
universe thermalizing with different values of will be characterized by different amplitudes 
of the cosmological density fluctuations. Examples of models of this sort have been given in 
Refs. 1,0. 

The application of the principle of mediocrity in our case will require comparing the 



^We use the term a priori in the sense that this distribution is independent of the existence of 
observers. 



expected numbers of civilizations in parts of the universe with different values of A and Q, 
which will be treated as random variables. In fact, for our purposes, it will be convenient 
to deal with an additional random variable, to- This is because one of the questions we are 
addressing is the coincidence (2), and galaxy formation can itself be modeled as a random 
process which takes place over a range of times for given Q and A. Instead of Q, it will be 
more convenient to use the density contrast on the galactic scale at the time of recombination, 
(Tree- Throughout the paper we assume that the universe is fiat, ^a + Qm = 1- 

The paper is organized as follows. We shall first consider the situation in which only the 
cosmological constant is allowed to vary, with all other parameters being fixed. In Section 
II we will show that the most likely values of A and to in this case are such that t\ ~ tfj. In 
Section III we shall argue that the most likely epoch for the existence of intelligent observers 
is to ~ tc- This completes the argument that coincidences (^ and (|]) are indeed to be 
expected in this class of models. In Section IV we discuss models where both A and arec 
are variable and outline the calculation of the probability distribution for tA and to- In 
our analysis of these models we go beyond the issue of the cosmic time coincidence and 
discuss the values of ^a and of the density contrast a^ec detected by typical observers. Our 
conclusions are summarized in Section V. 

II. WHY IS tA ~ tc? 

In this and the following section we assume that the cosmological constant A is the 



only variable parameter. Weinberg [|18[ was the first to point out that not all values of A 
are consistent with the existence of conscious observers. In a spatially fiat universe with 
a cosmological constant, gravitational clustering effectively stops at a redshift (1 + za) ~ 
(Pa/ PmoY^^, when pA becomes comparable to the matter density pm- (Here, pmo is the 
present matter density.) At later times, the vacuum energy dominates and the universe 
enters a de Sitter stage of exponential expansion. An anthropic bound on pA can be obtained 
by requiring that it does not dominate before the redshift Zmax when the earliest galaxies 



are formed 



PA<{1 + Zmax) PMO- (5) 

Weinberg took Zmax ~ 4, which gives pA ~ lOOpMo- 

One expects that the a priori probabihty distribution 'P*(pa) should vary on some char- 
acteristic particle physics scale, ApA ~ v'^- The energy scale rj could be the Planck scale 
rjpi ~ 10^^ GeV, the grand unification scale tjgut ~ 10^^ GeV, or the electroweak scale 
Vew ~ 10^ GeV. For any reasonable choices of rj and z^ax, ^Pa exceeds the anthropically 
allowed range of pA by many orders of magnitude. We can therefore set 

'P*(Pa) = const (6) 



in the range of interest |]T8[. With this fiat distribution, a value of pa picked randomly from 



an interval |pa| ^ p™ is likely to be comparable to p^ (the probability of picking a much 
smaller value is small). In this sense, the fiat distribution (|^) favors larger values of pa. 

The anthropic bound (^ specifies the value of Pa which makes galaxy formation barely 
possible. However, the principle of mediocrity suggests that we are most likely to observe 
not these marginal values, but rather the ones that maximize the number of galaxies. This 
suggests that A-domination should not occur before a substantial fraction of matter has 
collapsed into galaxies. The largest values of A consistent with this requirement are such 
that t\ ~ to- Hence, the coincidence (H) is to be expected if we are typical observers ||T9| . 

Let us now try to make this more quantitative. It will be convenient to introduce a 
variable 

x=S^ = sinh^ (-) , (7) 



where for convenience, we have defined ^a as the time at which ^a = sinh (1)^2^/ ~ 1.38^^/. 
At the time of recombination, for values of pA within the anthropic range, Xrec ~ Pk/ Prec ^ 1, 
where the matter density at recombination, prec-, is independent of A. We can therefore 
express the probability distribution for pA as a distribution for Xyec-, 
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ill yX'i'f^c) OC t'xXf^cjClX^^C') \ ) 

where v{xrec) is the number of galaxies formed per unit volume in regions with a given value 
of Xrec- The calculation of the distribution (|]) was discussed in detail by Martel et. al. 0. 
A simplified version of their analysis is given in the Appendix. 

Galaxies form at the time when the density contrast (evolved according to the linear 
theory) exceeds a certain critical value Ac(x). For small values of x, when the cosmological 
constant is negligible, we have Ac(x) ~ 1.69 as in the Einstein de Sitter model. However, 
it is known that Ac is slightly dependent on x, with Ac(oo) k, 1.63. Thus, Ac varies by 
no more than 4% in the whole relevant range, and in what follows we shall ignore its x- 
dependence. The number of galaxies wich have assembled up to a given time t for a given 
value of the cosmological constant (that is, up to a given x for a given value of x^ec) can 
thus be estimated as [^ 



v{x,Xrec) = erfc —1= . (9) 

\\/2arecG{x,Xrec)) 

The factor (^(x, x^ec) = ^reV^-^(^)) "where 

5 /l + xy/2 /-^ duj 

Q\ir) Jo a;i/6(i + ^)3/2' (10) 

accounts for the growth of the dispersion in the density contrast a on the galactic scale from 
its value a^ec at the time of recombination until time t{x). For small x we have F ^ x^^^, 
and perturbations grow as in the Einstein de Sitter model. However, at large x the growth 
of perturbations is stalled and we have F{oo) = (5/6)/3(2/3, 5/6) ~ 1.44. The number of 
galaxies that will assemble in a given interval of x will thus be given by 



dv{x^Xrec) OC CXp 



1 / Ac xlLl 



2 \F{x) ar, 



1 



rec 



F'jx) xli^ 



dx. (11) 



Multiplying by a fiat a priori distribution for Xrec, we have 

dV{x,Xrec) OC du{x, Xrec)dXrec- (12) 



The probability for an observer to live in a galaxy that formed in a given logarithmic interval 
of tc/tx can now be obtained by integrating (p!2| ) with respect to Xrec while keeping x fixed. 
The result is 

dV{tGlU) oc ol^F'F' f M n{talt^). (13) 

This distribution is shown in Fig. 1 (curve a). It has a broad peak which almost vanishes 
outside of the range .1 ^ (tG/t\) ^10. The maximum of the distribution is at tc/tA ~ 1.7 
and the median value is at tc/tA ~ 1.5. Thus, most observers will find that their galaxies 
formed at t ~ tA, and therefore the coincidence 

to ~ ^A (14) 

is explained. 

It is also of some interest to consider the distribution (|12]) without performing any inte- 
grations. By changing from the variables x and Xrec to the variables tc and ^a we have 



dV oc a:;^^ exp 



il(:i)"(!f)-(^)--«--^' (-) 



2F2 
where x = xitc/tA) and 

U = iAJarecf/\ec (16) 

is the time at which the density contrast on galactic scales would reach the critical value Ac 
in an Einstein-de Sitter model. Here we are not allowing for variations of dree, and therefore 
this time is just a constant. The probability density (^5|) per unit area in the (logtG,logtA) 
plane is plotted in Fig. 2. Note that the peak is in the region where t^ ~ ^a ~ ^o-- Different 
projections of this plot are useful. If we integrate along the vertical axis, then we obtain the 
probability distribution for the time when A dominates, which is equivalent to (|^), whereas 
if we integrate diagonally along (^g/^a) = const lines, we obtain (p^). 



dP/d[10g(tG/t^ )] 
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FIG. 1. The probability density per unit logarithmic interval oitc/tA, Eg. (|l3|) , is shown (curve 
a). The maximum is at tc/tA ~ 1.7 whereas the median value is at t^/^A ^ 1-5. The same 
distribution taking into account the cooling boundary tcb discussed in Section IV is also plotted 



(curve b). The parameters have been chosen so that tc6 = -5 t^ [see Eq. (32)]. 
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FIG. 2. The joint probability density ( [l5| ) per unit area in the plane log (tA/to-) (horizontal 

axis) log(tG/to-) (vertical axis), where ta is defined in Eq. ([l6|). 



III. WHY NOW? 



As we noted in the Introduction, one of the most puzzling aspects of the value of VL\ 
is related to the fact that the coincidence to ~ ih appears to be implying that we live in a 
special time. A similar problem exists even if a quintessence component [^ is assumed (see 
Section V). As we have shown in Section II, the epoch when giant galaxies are assembled, 
tc, is expected to roughly coincide with the epoch of cosmological constant dominance, ^a- 
Therefore, if we could explain why to ~ ^g, the puzzle of the cosmic age coincidence would 
be resolved. 

Most of carbon-based life may be expected to have appeared in the universe around the 
peak in the universal carbon production rate, at tcarbon- The main contributors to carbon 
in the interstellar medium are stars in the mass range 1-2 Mq, through carbon stars and 
planetary nebulae [^. Consequently, detailed simulations ^3[ show that the peak in the 
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cosmic carbon production rate is delayed only by less than a billion years compared to the 
peak in the cosmic star formation rate, tsFR, namely, 

if^carbon ~ ^SFR • (17) 

The appearance of intelligent life is further delayed by no more than a fraction of the 
main sequence lifetime of stars in the spectral range mid-F to mid-K ( 5-20 Gyr; [^,^). 
Following the main sequence, the expansion and increase in luminosity of stars spells the 
end of the possible existence of a biosphere on planets. Only stars in the above spectral 
range are expected to have continuously habitable zones around them (namely, ensuring the 
presence of liquid water and the absence of catastrophic cooling by CO2 clouds on planetary 
surfaces; [|^). Thus we have 

^IL ~ ^carbon ~ i^SFR • (18) 

The "present time" to can be defined as the time when a civilization evolves to the 
point where it is capable of measuring the cosmological constant and becomes aware of the 
coincidence (|I|). The experience of our own civilization suggests that, on a cosmological 
timescale, this time is not much different from tiL, 

to ~ tiL. (19) 

Carter [|r^ and others ||T^JT3[] used the principle of mediocrity to argue that the lifetime 
of our civilization is unlikely to be much longer than the time it has already existed, that is, 
~ 10^ yrs. If we are typical, then this should be the characteristic lifetime of a civilization. 
This would imply that Eq. ([19D is valid even if to is understood as the time when any 
astronomical observations can be made. Carter's argument has some force, but it is based 
on a single data point, and one may be reluctant to accept it, considering in particular its 
pessimistic implications. We note, however, that with our definition of to, Eq. (|T9|) is likely 
to be valid regardless of the validity of Carter's argument (that is, even if civilizations are 
likely to survive much longer than tn,). Combining (^) with (|IB|), we have that for a typical 
civilization 
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^0 ~ tsFR- (20) 

Finally, models of galaxy formation in hierarchical clustering theories propose that 
Lyman-break galaxies (at 2; ~ 3) are the first objects of galactic size which experience 
vigorous star formation [^. These objects therefore signal the onset of the epoch of galaxy 
formation, with cosmic star formation and galaxy formation being closely linked. In fact, 
the mergers and collisions of "sub-galactic" objects to produce galactic-size structures, are 
responsible for the enhanced star formation. In hierarchical models therefore 

to ~ tsFR ■ (21) 

The above relation is also supported by observations of the star formation history, show- 
ing that the star formation rate rises from the present to about 2; ~ 1, with a broad peak (of 
roughly constant star formation rate) in the redshift range z ~ 1-3 [^. This corresponds 
roughly to tsFR ~ ^o/3-to/8, in agreement with Eq. (^). In fact, more than 80% of the 
stars have already formed {Qgas/^stars ~ 0.18 p9|). 

Combining Eqs. (p!^) , (p0|) , (pT]) above we obtain the desired relation 

tp ~ tc ~ U . (22) 



IV. MODELS WITH VARIABLE A AND aREC 

In the previous discussion we have assumed a fixed value of the density contrast at 
recombination arec (or equivalently, a fixed value of Q). This determines the parameter 
to- = i^c/o'recY^'^trec appearing in the distribution (|T5|) and therefore, as it is clear from Fig. 
2, the most probable time at which the cosmological constant will dominate tA ~ ^o- 0- 

If CTree IS itsclf treated as a random variable, with a priori distribution V^{crrec)d^'^o'rec 
then the most probable value of ^a will of course have some dependence on "P*. However, as 
we shall argue, this dependence is not too strong provided that P* satisfies some qualitative 
requirements, in which case the most probable values of ^a and arec are actually determined 
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by the fundamental constants involved in the cooling processes which take place in collapsing 
gas clouds. 

A. The cooling boundary 

So far, we have assumed that all the galactic-size objects collapsing at any time form 
luminous galaxies. However, galaxies forming at later times will have a lower density and 
shallower potential wells. They are thus vulnerable to losing all their gas due to supernova 
explosions [|l^. Moreover, a collapsing cloud fragments into stars only if the cooling timescale 
of the cloud Tcooi is smaller than the collapse timescale Tgrav Otherwise, the cloud stabilizes 
into a pressure supported configuration |]35| , p!0| . The cooling rate of such pressure supported 
clouds is exceedingly low, and it is possible that star formation in the relevant mass range 
will be suppressed in these clouds even when they eventually cool. Hence, it is conceivable 
that galaxies that fail to cool during the initial collapse give a negligible contribution to z/. 
Fragmentation of a cloud into stars will be suppressed after a certain critical time which we 
shall refer to as the "cooling boundary" tcb |10l . 

To determine tcb, let us first consider the case of a matter dominated universe (not nec- 
essarily flat) without a cosmological constant. An overdensity which is destined to collapse 
can be described in the spherical model as a part of a closed Friedmann-Robertson- Walker 
(FRW) universe. The size of this spherical region at the time of recombination is such that 
it basically contains the mass of a galaxy. The virialization temperature and the density 
after virialization will be quite independent of what happens outside the region, depending 
only on its gravitational energy at the time tyir when it collapses. The virial velocity will 
then be given by Vmr ~ (GMg/LY^'^, where L is the size of the collapsing object at tmr- The 
density of the virialized collapsing cloud pmr is given by p!0| , |36[| 



p„. ~ W\Gtl,,)-\ (23) 

The virialization temperature can be estimated as T^ir ~ T^pvli^ ~ f^piG^ PvirMgY^^ . Here 
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Trip is the proton mass. The later an object collapses, the colder and more dilute it would 
be. 

If there is a cosmological constant, then these estimates still hold to good approxima- 
tion. Indeed, a spherical region will only collapse if its intrinsic "curvature" term is always 
dominant with respect to the cosmological constant term. The "potential" energy at the 
time of collapse and the properties of the virialized cloud will basically remain unaltered. In 
principle, a spherical region with a cosmological constant could enter a "quasistatic" phase 
where the gravitational pull is nearly balanced by the repulsion due to the cosmological con- 
stant. After a long period of time, this region might finally collapse and virialize to a large 
enough temperature. However, since the quasistatic phase is unstable we shall disregard 
this marginal possibility. 

The cooling rate t~Ii of a gas cloud of fixed mass depends only on its density and 
temperature, but as shown above both of these quantities are determined by t^r |0- The 
timescale needed for gravitational collapse is Tgrav ~ t^ij.. Therefore, the condition Tcooi < 
Tgrav givcs an uppcr bound tct on the time at which collapse occurs. Various cooling processes 
such as Bremsstrahlung and line cooling in neutral hydrogen and helium were considered in 
Ref. [|TU|. For a cloud of mass Mg ^ 10^^ Mq, cooling turns out to be efficient [^ for 



t < U ^ 3 • 10 V- (24) 

In any case, this value of tcb should be taken only as indicative, since the present status of 
the theory of star formation does not allow for very precise estimates. 

B. Likely values of ^a 

Let us now consider the probability distribution for the three independent variables 
X, Xrec and (Tree- This will be proportional to the number of galaxies forming at a time 
charachterized by x in a region with given values of arec and Xrec, 

pi ^1/3 



dl \X^ ^rec: ^rec) ^ ' *\^rec) G^P 



1 /A x^/'^V 



^ ^rec 
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"rec 



-^ ^rec 



dx dXrec dlnarec- (25) 



Let us assume for simplicity a power-law a priori distribution, 

V^{arec) oc o-^g", (26) 

where a is a constant. Then we can immediately integrate over arec and obtain 

Now we can integrate with respect to the "time" x at which galaxies assemble, from the 
time of recombination Xrec to the cooling boundary 

Xeb = sinh2(t,fe/U). (28) 

The integral is simply the difference in F° between the two boundaries in the integration 
range, and we shall neglect the contribution at Xrec- Finally, using ^a = trecXrel"^ "we obtain 
a probability distribution for ^a 

dV{t^) oc F°(sinh2(t,,/tA))tf "'rfln(tAAcfe). (29) 

Thus, the most probable value of ^a is determined by tcb and a. 

In Fig. 3, this distribution is plotted for different values of a ranging from 4 to 15. In 
all these cases we have 

tA ~ U- (30) 

The behaviour of the distribution is different for a < 3. Note that F{y) oc y^^'^ for small 
y, whereas F saturates at a constant value for large y. This means that if a < 3, the 
distribution (|2^) would favour very small values of ^a- The reason is that for a small a the 
a priori distribution is not too suppressed at large Urec, and it pays off to increase arec in 
order to obtain a large number of collapsed objects very soon after recombination. Therefore 
the time of A domination can be very short without interfering with galaxy formation. Of 
course, this would result in an overwhelming majority of the galaxies in regions which do 
not look anything like ours. On the other hand, if a > 3, small values of a^ec are preferred. 
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However, the value of arec should at least be large enough for galaxy formation to occur 
marginally before the cooling boundary tcb- Therefore, if the cosmological constant is not 
to interfere with galaxy formation, the result ( pOD is expected. More generally, we expect 
the relation (^) to be valid if the a priori distribution decreases faster than o"~^ at small 
arec- With tcb from Eq. ( pi[ ) and the ^a suggested by observations, the relation ( ppj ) is indeed 
satisfied. 



dP/d[log(t./t,J] 
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FIG. 3. The probability density per unit logarithmic interval of tj^/tf-b-, Eq. (E^), for different 
values of the parameter a (a = 4, 5, 10, and 15). 

C. Likely values of arec 

A probability distribution for arec can be obtained by integrating (|25|) first over Xrec over 
the relevant range sinh(a;.J/^tcfe/^rec) > a;^/^, and then over x. The result can be expressed as 



dV{(3) oc (3^^''-'^^'^G{(3) d\np, 



(31) 



where we have introduced the arec dependent parameter 
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tcb 



A \^/2t 



(Jr- 



tcb 



(32) 



and the function 

Gi/3) 



exp 



"2F2 Ua i 



4/3- 



F' 



2 Ua/ 



4/3- 



dx. 



(33) 



The function G'(/3) is plotted in Fig. 4 (thick sohd hne). It stays constant for /3 < 1 (towards 
large arec) and it drops to zero around P ~ 10. For larger /? it falls off as /?^/^exp(— /5^/^/2). 
This function is multiplied in (|3TD by the factor o"^~", which depends on the a priori distri- 
bution for a^ec- If this factor is a decreasing function of arec (i-e. a > 3) then (^TJ) peaks 
between 1 < jS ^ 10. This is illustrated in Fig. 4 (thin curves) for a = 4 and 5. From the 
definition of (3 we have 



A. 



2/3 



(J^ 



1.1-10-^/5 



-2/3 



(34) 



(1 + Zrec) \ 'iPHotcb V^M J 

Here, we have used the relation Hx^^^ = 2/(3tAV^M), where all quantities (including the 
matter density parameter VLm) are evaluated at the present time, and for our numerical 
estimate we have taken z^ec = 1000, tc6 = 3- lO^'^yr, Hq = 100 h km s~^ Mpc~^, with h = .7, 
and Qm = -3. For /3 ~ 1, as suggested by the distribution (^Tj), we have that the most likely 
values of arec are of the order of 10^'^. This is close to the observationally suggested values 
a,ec = (3.3 -2.4)10-3 g. 
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FIG. 4. The probability density per unit logarithmic interval of /? oc arec , Eq. (|3l[), for a = 4 

and 5 (thind solid lines). The function G(/?) is represented by a thick solid line. 

Anthropic bounds on the density contrast have been recently discussed by Tegmark and 
Rees [|l^. Instead of arec, they used the amplitude of the density fluctuations at horizon 
crossing, Q; the relation between the two is roughly Q ~ lO^^a^ec- They imposed a lower 
bound on Q by requiring that galaxies form prior to the "cooling boundary", t^- ^ tct- This 
gives Q ^ 10^^. To obtain an upper bound, it has been argued [pOi|TO| that for large values 
of Q galaxies would be too dense and frequent stellar encounters would disrupt planetary 
orbits. To estimate the rate of encounters, the relative stellar velocity was taken to be the 
virial velocity v^r ~ 200 km s~^, resulting in a bound Q ^ 10~^. However, Silk [^ has 



pointed out that the local velocity dispersion of stars in our Galaxy is an order of magnitude 
smaller than Vmr- This gives Q ^ 10~^, which is a rather weak constraint. This issue does 
not arise in the approach we take in the present paper, since in our case large values of Q 
are suppressed by the a priori distribution P*(crrec)- 



D. The time coincidence 

Finally, we should check that the introduction of a cooling boundary does not spoil 
the coincidence tc ~ ^a- In fact, this seems rather clear from Fig. 2. Introducing the 
cooling boundary basically amounts to disregarding the probability density above a certain 
horizontal line to = tab- The probability distribution for tG/t\ below the horizontal line 
is somewhat different from that in the whole plane, but clearly it still peaks around to ~ 
^A- To quantify this effect, we have integrated (|T^ with respect to Xrec over the range 
smh(xmtcb/trec) > x^l"^ . The resulting distribution for x is proportional to the integrand in 
the right hand side of (|33|) . For /? = 2, this probability density is shown in Fig. 1 (curve 6). 
The peak is only slightly shifted towards smaller values of tcltk- 

Cooling failure is not the only mechanism that can in principle inhibit the number of 
civilizations at low o^ec- It is possible, for example, that the stellar initial mass function 
(IMF) depends on the protogalactic density pmr-, so that the number of carbon forming 
stars drops rapidly towards very low values of pmr- If the a priori distribution V^iorec) 
is a decreasing function of Orec-, this can result in a peaked distribution dV /d\n.t\. Quite 
similarly, if the number of relevant stars grows towards smaller p^r, a peaked distribution is 
obtained for an increasing function ^^(crrec)- Our present understanding of star formation 
is insufficient to determine the dependence of the IMF on p^j^, but once it is understood, 



the probability distribution for ^a can be calculated as outlined above |^ . 



V. CONCLUSIONS 

In this paper we suggested a possible explanation for the near-coincidence of the three 
cosmological timescales: the time of galaxy formation tc-, the time when the cosmological 
constant starts to dominate the energy density of the universe ^a, and the present age of the 
universe to- Since this coincidence involves specifically the time of our existence as observers, 
it lends itself most naturally to the consideration of anthropic selection effects. 
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We considered a model in which the cosmological constant is a random variable with a 
flat a priori probability distribution. We showed that a typical galaxy in this model forms 
at a time to ~ ^a- We further demonstrated that a typical civilization should determine 
the value of the cosmological constant at to ~ ^g- Thus we should not be surprised to find 
ourselves discussing the cosmic time coincidence. 

We also considered a model in which both the cosmological constant A and the density 
contrast a-rec are random variables. The galaxy formation in this case is spread over a much 
wider time interval, and we had to account for the fact that the cooling of protogalactic 
clouds collapsing at very late times is too slow to allow for efficient fragmentation and star 
formation. We therefore disregarded all galaxies formed after the "cooling boundary" time 
tcb- We assumed that the a priori distribution for arec is a decresing power law and found 
that, for a sufficiently steep power, a typical observer detects arec ~ 10^^ — 10~^, close to 
the values inferred from observations. Such observers are likely to find themselves living at 
to ~ tA in a galaxy formed at t^ ~ tA in a region of the universe where tA ~ tcb, also close 
to the observationally suggested value. 

Our model with variable A and arec can be developed further in several directions. Instead 
of taking a flat distribution for pA and a power-law distribution for arec, one could use the 
methods of Refs. ||T7| , ^ to calculate the a priori distributions for these variables in the 
framework of some inflationary model. One could also use a more reflned model of structure 
formation and improve on our treatment of cooling failure, replacing the sharp cutoff at 
t = tcb with a more realistic model. We believe, however, that even in the present, simplifled 
form our model indicates that an anthropic selection for A and arec is a viable possibility. 

Finally, we should note that the coincidence in the timescales requires an explanation 
even in models involving a quintessence component [^. In models of quintessence the 
universe at late times is dominated by a scalar fleld 0, slowly evolving down its potential 
V{(f)). It has been argued (by Zlatev et al. [0) that such models do not suffer from the 
cosmic time coincidence problem, because the time t^ of 0-domination is not sensitive to 
the initial conditions. This time, however, does depend on the details of the potential V^(0), 
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and observers should be surprised to find themselves living at the epoch when quintessence 
is about to dominate. More satisfactory would be a model in which the potential depends 
on two fields, say (f) and Xi with x slowly varying in space, making the time of ^-domination 
position dependent. Such models are not difficult to construct in the context of inflationary 
cosmology. One could then apply the principle of mediocrity to determine the most likely 
value of t^. 
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APPENDIX: THE PROBABILITY DISTRIBUTION FOR A 

In this Appendix we briefly discuss the probability distribution for the cosmological 
constant, giving a simplified version of the calculation presented in 0. 

In a universe where the cosmological constant is non- vanishing, a primordial overdensity 
will eventually collapse provided that its value at the time of recombination exceeds a certain 
critical value 6^^'^. In the spherical collapse model this is estimated as 5'^^^^ = 1.13 xHl (see 
e.g. [0). Hence, the fraction of matter that eventually clusters in galaxies can be roughly 



approximated as pO| , p! 



Here, erfc is the complementary error function and arec{Mg) is the dispersion in the density 
contrast at the time of recombination on the relevant galactic mass scale Mg ~ IO^^Mq. 
The logarithmic distribution dV /din Xrec = Xreci^{xrec) is plotted in Fig. 5. It has a rather 
broad peak which spans two orders of magnitude in Xrec, with a maximum at 

a;^:f ^ 2.45 al,. (36) 
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In accordance with the principle of mediocrity, we should expect to measure a value of the 
cosmological constant within this broad peak of the distribution. And indeed, this may 
actually be the case. The distribution (|35| ) is characterized by the parameter area- As 
noted by Martel et al. 0, this parameter can be inferred from observations of the cosmic 
microwave background anisotropics, although its value depends on the assumed value of the 
cosmological constant today. For instance, assuming that the present cosmological constant 
is f2A,o = -8, and the relevant galactic co-moving scale is in the range R = {l — 2)Mpc, Martel 
et al. found a-rec = (2.3 — 1.7) 10^'^. In this estimate, they also assumed a scale invariant 
spectrum of density perturbations, a value of 70km s~^ Mpc~^ for the present Hubble rate, 
and they defined recombination to be at redshift Zrec ~ 1000 (this definition is conventional, 
since the probability distribution for the cosmological constant does not depend on the 
choice of reference time). Thus, taking into account that x scales like (1 + z)~^ in equation 
(^), one finds that the peak of the distribution for the cosmological constant today is at 
Xg^" ^ 29.8 — 12. The value corresponding to the assumed ^a^ = .8 is Xq = 4, certainly 
within the broad peak of the distribution and not far from its maximum. If instead we assume 
that the measured value is Qa,o = -7, which corresponds to Xq = 2.33, the new inferred values 
for (Tree corrcspoud to the peak value Xq^"" ~ (88—34). In this case, the measured value would 
be at the outskirts of the broad peak, where the logarithmic probability density is about an 
order of magnitude smaller than at the peak, but still significant. Thus, even though there 
may be uncertainties in the inferred value of arec on the relevant scales, it seems fair to say 
that any observed value of ^a.o ^ -7 is in good agreement with the principle of mediocrity. 
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FIG. 5. The probability density per unit logarithmic interval of Xrec^^rec- 
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